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In supersymmetric models neutrino data can be explained by i?-parity violating operators 
which violate lepton number by one unit. The so called bilinear model can account for the 
observed neutrino data and predicts at the same time several decay properties of the lightest 
supersymmetric particle. In this paper we discuss the expected precision to determine these 
parameters by combining neutrino and LHC data and discuss the most important observables. 
We show that one can expect a rather accurate determination of the underlying i?-parity 
parameters assuming mSUGRA relations between the i?-parity conserving ones and discuss 
briefly also the general MSSM as well as the expected accuracies in case of a prospective e+e~ 
linear collider. An important observation is that several parameters can only be determined 
up to relative signs or more generally relative phases. 
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^ I. INTRODUCTION 

(N 

,— , With the start of the LHC the exploration of the terascale has begun and, thus, in the search 

^ for extensions of the Standard Model (SM) significant higher mass scales can be tested as in the 

^ past. Supersymmetry (SUSY) is among the most favoured candidates as it allows for example the 

unification of gauge couplings and stabilizes the hierarchy of the Planck scale (or the GUT scale) 
, and the electroweak scale. Once supersymmetry has been found, there will be two main tasks: (i) 

to determine the underlying parameters of the model and (ii) to check if the minimal model is 
^ realized in nature or an extended one. 

OO In principle one can supersymmetrize all known mechanisms to generate neutrino masses and 

^2 mixing angles, see e.g. ref. [ ]. However, supersymmetry offers an intrinsic possibility to generate 

neutrino masses, namely i?-parity violation (RPV), for an review see e.g. [ ]. In this case neutrino 
masses are generated either via mixing with neutralinos or via loop effects [ - ]. We will be 
2 constraining ourself to bilinear R-parity violation (BRpV) [ - ] , thus breaking lepton number but 

^-H not baryon number. This model can be viewed as the effective theory of a spontaneously broken 

^ i?-parity [9-12]. The BRpV neutrino mass model has only a few free parameters and therefore is 

very predictive. Furthermore, in contrast to trilinear RPV neutrino mass models, the constraints 
^ from LEP on the i?-parity violating couplings are automatically satisfied as the couplings are small 

due to the requirement of explaining correctly neutrino data. 

The RPV couplings giving rise to neutrino masses are also responsible for the decay properties of 
the neutralino. Therefore, an important smoking gun signal for these models is the strong connection 
between neutralino physics and neutrino mixing parameters. Some ratios of the branching ratios 
of the neutralino, or more generally the lightest supersymmetric particle, are related to neutrino 
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mixing angles [13-16]. In particular, approximately the same number of muons as taus are expected 
along with a VF-boson because their ratio is given by tangent of the nearly maximal atmospheric 
mixing angle. By measuring the decay properties of the neutralino, which is likely to be done by 
the LHC, a severe test of this model is possible. 

Due to the smallness of the RPV couplings the neutralino will have a long lifetime, but short 
enough for it to mainly decay within the detectors at the LHC. The prospects for collider discovery 
of RPV, responsible for the neutrino masses and mixings, in mSUGRA have been thoroughly 
studied [ - ]. 

Once this scenario is confirmed an important question will be how well the underlying parameters 
can be determined. There have been several studies within the MSSM with conserved i?-parity 
[22-28]. In this paper we want to focus on the determination of the RPV parameters taking into 
account neutrino and LHC data. We will focus on a specific scenario for the details but comment 
on how this can be extrapolated to others. We will also work out which are the most sensitive 
observables. 

This paper is organized as follows: in the next section we will recall briefly the main features 
of the BRpV model. In Section HI we discuss the details of the fit procedure and present the fit 
results in Section IV and then draw in Section V our conclusions. 



II. EXPLICIT BILINEAR i?-PARITY VIOLATION 

The MSSM with explicit bilinear i?-parity violation is specified by the superpotential ["^] 

WbrpV = Wmssm + eiZiHu , (1) 

where the last term explicitly violates both i?-parity and lepton number in all three generations. In 
addition one has to add corresponding terms to the soft SUSY breaking potential, 

T .BRpV _ T^MSSM \ TD ^ J TT 

which induce vacuum expectation values (VEVs) Vi (z>^) for the sneutrinos. Using the tadpole 
equations one can take the sneutrino VEVs as input instead of the Bi []. 

One important aspect of this model is that the i?-parity breaking terms give rise to mixings 
between SM and SUSY particles. In the neutral fermion sector the mixing between neutralinos and 
neutrinos leads to one massive neutrino at tree-level while the other two neutrinos acquire masses 
through loop corrections [2, 7, b]. The effective neutrino mass matrix at tree- and one-loop level is 
given by 

(^efr,Lo)ij = aKiKj (3) 
(^eff,NLo)ij = hkikj + c{Kiej + eikj) + deiej , (4) 

where a, 6, c, d are functions of i?-parity conserving parameters and are the so called alignment 
parameters 

Ai = iJiVi + Vdti . (5) 

For example a is given by 

-{g^Mi + g'^M2) 



(6) 
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and b is equal to a plus radiative corrections. From the mixing matrix that diagonalizes eq. (4) one 
obtains expressions for the neutrino mixing angles in terms of the i?-parity breaking parameters 
which can approximately be expressed as 



tan 6>23 - ' ^ ^2 _^^2 ' - (^^ j , (7) 

where = Vtree,zj^j and Vtree diagonalizes the tree-level neutrino mass matrix eq. (3). As one can 
see from eqs. (4) and (7) both the neutrino masses and the mixing angles are predicted in terms 
of the i?-parity breaking parameters. It turns out that the approximation for tan^ 023 is relative 
insensitive when going from tree-level to one-loop level. The one for tan^ Qy^ can be quite sensitive 
to loop corrections, in particular if e2A2e3A3 > 0. [7, 8]. This will also manifest itself later in the 
fits discussed in Section IV. The most important loop contributions are due to sbottom-bottom 
and stau-tau loops [ , ] which are in both cases proportional to eje/e/|/ip. However, there are also 
regions in parameter space where the mixing between sleptons and the charged Higgs boson [8] 
and/or the sneutrino-neutrino loops give large contributions [ ]. 

The violation of i?-parity also implies that the lightest supersymmetric particle (LSP), which we 
assume to be the lightest neutralino here, is no longer stable and typically decays inside the detector 
once the parameters are adjusted to satisfy the neutrino constraints [14-16, 31]. The parameters 
that determine the decay properties of the LSP are the same parameters that lead to neutrino 
masses and oscillation which implies that there are correlations between the neutralino branching 
ratios and the neutrino mixing angles [lo, i^], e.g. 

This can be most easily seen by performing first an approximate diagonalization of the chargino and 
neutralino mass matrices as in [32] . The parts of the corresponding mixing matrices responsible 
for the mixing between neutrinos and neutralinos as well as charged leptons and charginos can be 
expressed in terms of 



A7 A,- 

(9) 



and enter also the corresponding couplings in the decays used in eq. (8). Moreover, the decay length 
of the LSP is proportional to the i?-parity breaking parameters [i4, 19, 21]. All these facts can be 
used to determine these parameters once information on the i?-parity conserving ones is available. 



III. FIT PROCEDURE 

LHC will provide first information on the SUSY parameters once the corresponding signals are 
observed. However, it will be unlikely that the complete spectrum will be discovered and, thus, 
the first parameter fits will be performed within specific high scale models [ - ]. Therefore, the 
fits presented in this paper are performed in mSUGRA^ models which are augmented by bilinear 



^ Taking mSUGRA is not crucial, as the explanation of neutrino data does not depend on this assumption and can 
equally well be explained in GMSB [31], AMSB [33-35] or the general MSSM [. _]. 
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i?-parity breaking parameters at the electroweak scale. ^ These models therefore have eleven free 
parameters, namely the five mSUGRA parameters mo, tan/3, and sgn(/x) and the six 

bilinear i?-parity breaking parameters e^, and (or Vi respectively). For the theoretical predictions 
of the neutrino oscillation data, the LSP decay properties, and LHC/ILC observables, SPheno^[36] 
version 3.0.beta54 has been used. 

In order to measure the agreement between the data and the model for a particular choice of 
parameters a simple function is used, 

x'(a) = E(yi-/iH)V^f, (10) 

i 

where the i/i are data points with their associated uncertainties or experimental errors cr^ and the 
fi{a) are theoretical predictions for these data points at the point a in parameter space. The data 
points used for the fits were also calculated by SPheno for a specific mSUGRA point, where the 
i?-parity breaking parameters are not explicitly specified but are calculated iteratively such that 
the predicted neutrino mixing angles and squared mass differences lie in the 3a confidence region as 
given in [37, Table Al]. If the i?-parity breaking parameters as determined by SPheno are denoted 
by a, then the data points are equal to their predictions at the point a where the i?-parity breaking 
parameters are explicitly specified, i.e. i/i = 

Since the data points used in the fits are themselves theoretical predictions, the absolute minimum 
of eq. (10) is trivially given by Xmin X^(S) — 0- Instead of finding the parameter point with the 
best goodness of fit, the purpose of data fitting is then to estimate the parameter errors based on 
the uncertainties cr^, or to locate other minima in a multimodal x^-l^ndscape that have equally good 
X^- values. In order to find all minima in a specific region of parameter space and to determine the 
boundaries of x^(^) ^^e error estimation, the Minuit Migrad[ ] optimization algorithm was 
used repeatedly at random starting points in parameter space. This procedure gives a good coverage 
of the x^-lctndscape and finds with high probability all minima which lie in the region that is 
bounded by the starting points. As interface between SPheno and Minuit the general purpose fitting 
program Kaimini^, which provides different deterministic and stochastic optimization algorithms 
and works with any program that implements the SUSY Les Houches Accord (SLHA)[39, 40], is 
used. 

IV. RESULTS 
A. Fit Setup 

Various fits with different free parameters and data points were carried out. The following 
subsections discusses fits where the i?-parity breaking parameters are free parameters and different 
observables that depend on these parameters (that are the neutrino oscillation data and the Xi 
decay properties) are used as data points. For each of these combinations three fits are performed 
which differ with respect to the mSUGRA parameters being fixed or free parameters of the fit 
and the set of corresponding data points. In the first setup the mSUGRA parameters are fixed 



^ This model is sometimes called RmSUGRA or BRpV-mSUGRA in the literature. 

^ The latest SPheno version can be obtained from: http://physik.uni-wuerzburg.de/~porod/SPheno.html 
^ The latest Kaimini version can be obtained from: https://github.com/fthomas/kaimini 
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and only the neutrino and/or neutralino observables are used as data points. In the second and 
third setup the mSUGRA parameters mo, and tan/3 are free parameters in addition to 

the i?-parity breaking parameters. The additional data points of the second setup are the "edge 
variables" (m2)«dge^ (m2„)«dge^ (m^ ,)thres, (ml^)*^-, (m^ )^dge, and (m^J^^dge^ ^^ere the predicted 

relative uncertainties were taken from [ , Table 5.13]. This setup is denoted as "LHC". The third 
setup includes in addition to the edge variables the masses of the Xii xt^ slept ons, and the ii 
assuming that are measured at a prospective future e+e~ linear collider such as ILC or CLIC. The 
relative uncertainties of these observables were taken from [ , Table 5.14]. This setup is denoted as 
"LHC+ILC". To be conservative the total uncertainties in both setups were obtained by summing 
statistic and systematic uncertainties linearly. We note, that taking the relative uncertainties for 
the collider observables equal for different points of parameter space is a strong assumption which 
would have to be confirmed by individual studies which are however not available in the literature. 
However, as we will see below the uncertainties on the RPV parameters are dominated by the 
uncertainties of the measurements of neutralino decay branching ratios and that the uncertainties 
on the measurements of masses and edge variables are only sub-dominant. 

Fits for the various setups were performed for different SUSY benchmark points, including 
SPS la^ [ ], la, lb, and 3 [ ]. In this paper we discuss in detail as a typical example SPS la' 
working out the main features as most experimental studies for the collider observables used have 
been performed for points close by and, thus, the assumption of using the same relative uncertainties 
can be more easily justified. The results for the other study points and additional plots for SPS la' 
are given on our web page [oc-x]. We are aware that SPS la' is potentially excluded by recent LHC 
data [45, 46] although the corresponding searches have been performed for the i?-parity conserving 
case. However, we have checked that qualitative features do not depend on the point under study, 
e. g. they are the same for SPS 3 which is not excluded by existing data. 

B. Neutrino oscillation data 

The free parameters of these fits are the six i?-parity breaking parameters while the data points 
are the neutrino oscillation data, which are the two mass squared differences Am^^^ and Amg^j and 
the three mixing angles tan^ ^atm, tan^ ^goi, and sin^ ^13. As uncertainties for the data points the 
mean value of the upper and lower la errors from Table Al of [ ] were used. To get meaningful 
results from the fits, the number of degrees of freedom, which is the number of data points minus 
the number of parameters, must be equal to or greater than zero. Therefore the Xi decay width was 
used as additional data point and we assumed that it can be measured with an accuracy of 15 % 
[47] . We perform this fit separately to get an understanding how the different sectors contribute to 
a global fit. 

One general result of these and all following fits is that x^(^) is symmetric under the transfor- 
mations —Ai and — e^, i. e. for a given minimum there are several other minima with the 
same goodness of fit which only differ in the signs of the i?-parity breaking parameters. Independent 
of the specific mSUGRA point, the x^-lctndscape of these fits is rather complex with 8 minima per 
sign combination of the A^, resulting in 64 minima in total and each of them having x^-values of less 
than 10~^. For a specific choice of signs for the A^, the 8 minima divide into two classes of 4 minima 
each which differ by the importance of their individual loop contributions. These minima can be 
differentiated by the sign of the ratio e2e3/(A2A3) [7, 8]. For e2e3/(A2A3) < 0, which is realized 
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Figure 1: Likelihood maps of fits at SPS la' for various parameter combinations using neutrino data and the 
neutrahno decay width: 61-62 (upper left), 6^-|A| (upper right), A1-A3 (lower left) and A2-A3 (lower right). 
These maps are projections of the six-dimensional parameter space onto the corresponding planes such that 
points with a higher likelihood cover points with a lower likelihood. The regions with maximum likelihood 
(which are colored black) correspond to minimal x^-values. See section IV B for further details. 

in the lower two plots of Figure 1 for A3 < 0.05 GeV^, the most important contributions to the 
solar masses and mixings come from bottom/sbottom and chargino/charged-scalar loops, while for 
minima where e2e3/(A2A3) > (regions with A3 > 0.05 GeV^) the corrections from bottom/sbottom 
loops are dominant. In the latter case also eqs. (7) get sizable corrections, in particular tan^0i3. 
This can also be inferred from the right upper plot of Figure 1. In the region in the upper left 
corner the atmospheric neutrino mass scale are essentially given by the tree level whereas in the 
lower two regions sizeable loop corrections are needed. The ratio e^/|A| gives the importance of 
the loop contributions relative to the tree-level-induced neutrino masses [ ] . In the left upper plot 
we show preferred regions in the 62-63 plane to demonstrate the sign ambiguities of the underlying 
parameters. Note, that in the lower two plots only one quadrant is shown and the remaining ones 
can be obtained by mirroring similar to the upper left. 

The absolute uncertainties on these parameters depend clearly on the i?-parity conserving 
parameters but the relative uncertainties are fairly insensitive to the i?-parity conserving parameters. 
However, up to now we have ignored their uncertainties as we only took neutrino data into account. 
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Figure 2: Likelihood map of the e2A2-e3A3 plain of fits at SPS la' with the Xi decay properties (compare 
subsection IV C). The plot shows that the distinct regions which correspond to different sets of open decay 
channels have different signs for €2^2 and 63A3 and that for the global minimum €2^2 < and e3A3 > 0. A 
fit with the neutrino and neutralino data combined will therefore favor points with this sign combinations. 



In the next sections we will add collider observables to allow also for a variation of these parameters. 



C. Neutralino decay properties 

We now investigate how the neutralino decay properties can be used to get information on the 
i2-parity breaking parameters without using neutrino data. We will use the Xi decay width and 
branching ratios as data points for the fits but keep the other collider observables fixed. Of the 
branching ratios only those were taken into account whose values exceeded 0.01 %. The relative 
uncertainties of the branching ratios were assumed to be ABR^/BRi = 2/y/N~i with TV^ = 2 • 10^ BR^ 
for a LHC integrated luminosity of 100 fb~^ [ ]. In the uncertainties of branching ratios whose 
final state contained quarks or r leptons we reduced the corresponding number Ni to A^i/10 in 
order to take potential uncertainties in the jet reconstruction into account. In principle this has to 
be checked by detailed Monte Carlo studies which are however beyond the scope of this article. 

As already mentioned in subsection IV B, the x^(^) of this fit is also symmetric under sign 
transformations of the A^, i. e. also the decay properties of the Xi do not depend on the sign of the 
alignment parameters. However, as can be seen in Figure 2 looking at the parameter combinations 
62^2 and eaAa one can identify a preferred quadrant, in this case the upper left, where all points 
with likelihood larger than 0.9 are located. 

Compared to the previous fits, the regions with a high goodness of fit is reduced which is partly 
due to the higher number of degrees of freedom of this fit (13 data points and 6 free parameters). 
There are now also distinct regions with a high goodness of fit (each with its own local minimum) 
that correspond to different sets of decay channels whose branching ratio exceeds the aforementioned 
threshold of 0.01 %. These regions can be seen in the two right plots of Figure 3. Within this setup 
the most important observables for the determination of the i?-parity breaking parameters can 
be inferred by counting the frequency of those observables with the highest x^-contributions for 
different minima. For SPS la^ one finds for minima with x^ < 1 that the four branching ratios with 
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Figure 3: Likelihood maps of fits at SPS la' for various parameter combinations using neutralino decay width 
and neutralino branching ratios: 61-62 (upper left), e^-|A| (upper right), A1-A3 (lower left) and A2-A3 (lower 
right). In the lower two only one quadrant is shown, the others can be obtained by mirroring with respect to 
the axes. These maps are projections of the six-dimensional parameter space onto the corresponding planes 
such that points with a higher likelihood cover points with a lower likelihood. The regions with maximum 
likelihood (black) correspond to minimal x^-values. See section IV C for further details. 

the highest x^-contributions are 

BR(x? ^rri^), BR(x? ^ rH , BR(x? ^ , BR(x? ^ &H . (11) 

Note, that this finding depends to some extent on the parameter point under study, e. g. for SPS lb 
the most important branching ratios are 

BR(x? ^ r/iz.) , BR(x?^rei/), BR(xS ^ 5^) , BR(xS ^ 6H , (12) 

whereas for SPS 3 they are 

BR(xS^rezy), BR(xS ^ /^H , BR(xS ^ , BR(xS ^ &H • (13) 



Moreover, it is obvious by comparing Figures 1 and 3 that neutrino data and branching ratios give 
complementary information. 
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Figure 4: Likelihood maps of fits at SPS la' for various parameter combinations using neutralino decay 
width and neutrahno branching ratios: 61-62 (upper left), 6^-|A| (upper right), A1-A3 (lower left) and A2-A3 
(lower right). In all but the upper right only one quadrant is shown, the others can be obtained by mirroring 
with respect to the axes. These maps are projections of the six-dimensional parameter space onto the 
corresponding planes such that points with a higher likelihood cover points with a lower likelihood. The 
regions with maximum likelihood (black) correspond to minimal x^-values. See section IV D for further 
details. 



D. Combination of neutrino and neutralino data 

In this section we are combining the previous two fits and add also other observables: the 
i?-parity breaking parameters are now fitted to the neutrino oscillation data and to the Xi decay 
properties. Due to the combination of the data points, the regions with a high goodness of fit is 
significantly reduced as can be seen in Figure 4. Note, the different scaling of the axes compared to 
the previous plots. For SPS la^ the parameter point a and all of its reflections with respect to the 
A^ axes are now the only minima with — 0. For these points €2^2 < and e3A3 > hold. In the 
range < 10 there is only one additional minimum with x^ ^ 1-5 and for which 62 A2 > and 
e3A3 < 0. 

Tables I-III give the uncertainties of the i?-parity breaking parameters inferred from these fits 
for the "neutrino+neutralino data", the "LHC", and the "LHC+ILC" setup for SPS la' and SPS 3. 
The main uncertainties are due to the uncertainties of the neutrino data and the neutralino decay 
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SPS la' 



SPS 3 



parameter 



best fit 



2(7 3(7 



ei [10-2 GeV] 
€2 [10-2 GeV] 
£3 [10-2 GeV] 



Ai 

A2 

As 



IQ-i GeV' 
10-2 GeV^ 



10-2 GeV' 



3.87 ±0.03 ±0.07 to\l 



-5.14 ±0.04 ±0.09 ±0.13 


^•^^-0.08 


+0.22 
-0.19 


+0.32 
-0.29 


03+0-04 

'-'•'-''^-0.10 


+0.08 
-0.13 


+0.11 
-0.16 


O.'±J-_0.04 


+0.11 
-0.10 


+0.16 
-0.15 


6.41 ±0.05 


+0.12 
-0.13 


±0.18 



parameter 



best fit 



2a 3a 



ei [10-2 GeV] 

62 [10-2 GeV] 

63 [10-2 GeV] 
10-1 GeV^ 



Ai 

A2 



10-2 GeV' 



Aa 10-2 GeV' 



6.88 ±0.04 toil ±0.18 

-9.14 ±0.05 ±0.17 ±0.22 

8.09 ±0.10 ±0.33 toil 

n n^+O-02 +0.05 +0.O6 

'-'•'-'^-0.11 -0.15 -0.17 



10.84 ±0.06 ±0.19 ±0.24 



10.84 ±0.06 ±0.21 toi? 



Table I: Parameter point a at SPS la' and SPS 3 with l<j, 2cr, and 3a uncertainties (12 d. f.) from fits with 
the neutrino oscihation data and the Xi decay properties. 



SPS la' 



parameter 



best fit 



2a 3a 



61 [10-2 GeV] 

62 [10-2 GeV] 

63 [10-2 GeV] 
10-1 GeV^ 



Ai 
A2 



3.87 ±0.03 t'ofr -oil 
-5.14 ±0.04 ±0.09 ±0.15 



4.55 



+0.10 

-0.08 



+0.22 +0.33 
-0.19 -0.30 



10-2 GeV' 



A3 10-2 GeV' 



n no+0.04 +0.07 +0.10 
'-'•'-''^-0.09 -0.14 -0.16 



6.41 ±0.04 ±0.10 ±0.16 



6.41 ±0.05 ±0.12 toil 



SPS 3 



parameter 



best fit 



2a 3a 



61 [10-2 GeV] 

62 [10-2 GeV] 

63 [10-2 GeV] 
Ai [lO-i GeV^ 
A2 [10-2 GeV^ 
A3 [10-2 GeV^ 



6.88 ±0.05 ±0.12 ±0.20 
-9.14 ±0.06 ±0.15 ±0.25 
8.09 ±0.11 ±0.26 toil 

n n^+O-02 +0.06 +0.09 
'-'•'-'^-0.12 -0.15 -0.20 

10.84 ±0.06 ±0.15 toit 



10.84 



+0.07 +0.15 +0.25 



-0.06 



-0.17 -0.26 



Table II: Parameter point a at SPS la' and SPS 3 with la, 2cr, and 3a uncertainties (30 d. f.) from fits with 
the "LHC" setup. 



branching ratios a s can bee seen by comparing Table I with Tables II and III. Except from Ai 
the percentage errors of the i?-parity breaking parameters at the 3a-level is less than 8 %. Ai is 
particularly difficult as the corresponding branching ratio Xi ^ VF^e+ is very small leading to the 
large uncertainties obtained. Note, that within the top-down approach, this means fitting the high 
scale mSUGRA parameters, there is hardly an effect on the uncertainties of the RPV parameters 
when going from the LHC to ILC. However, we do not expect significant differences when a fit of 
the i?-parity conserving parameters is performed at the electroweak scale provided the neutralino 
and chargino masses can be measured within 1 — 2 % which is feasible at the ILC [41, 48] or at 
CLIC [ ]. For the masses of the staus, sbottoms and the corresponding mixing angle one would 
need an accuracy of 10 — 20 % percent to keep the corresponding uncertainties in the calculation of 
the neutrino masses and mixing angles small enough so that they are sub-dominant compared to 
the uncertainties due to the branching ratios uncertainties. An important question will of course be 
on how much better the neutralino branching ratios can be measured at the ILC or CLIC compared 
to LHC. 
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parameter 



SPS la^ 
best fit 



2cr 



ei [10-2 GeV] 
€2 [10-2 GeV] 
63 [10-2 GeV] 
Ai [lO-i GeV^ 
As [l0-2 GeV^ 
10-2 GeV^ 



3a parameter 



+0.08 
-0.07 



3.87 ±0.03 
-5.14 ±0.04 ±0.09 



4.55 



+0.10 
-0.08 



os+°-°'^ 

'-'•^'^-0.09 



6 41+0-04 

'^•^-'--0.05 



6.41 ±0.05 



+0.22 
-0.19 

+0.07 
-0.13 



±0.10 



+0.10 
-0.11 



±0.12 
±0.15 

+0.34 
-0.31 

+0.10 
-0.16 

+0.17 
-0.16 

+0.18 
-0.17 



SPS 3 
best fit 



2cr 3cr 



ei [10-2 GeV] 
€2 [10-2 GeV] 
63 [10-2 GeV] 
10-1 GeV^ 
10-2 GeV^ 



Ai 

A2 



10-2 GeV 



6.88 ±0.05 ±0.12 ±0.20 
-9.14 ±0.06 ±0.15 toil 
8.09 ±0.11 ±0.26 toil 
n n^+o-02 +0.05 +0.09 

'-'•'-'^-0.12 -0.15 -0.19 

10.84 ±0.06 ±0.15 toil 



-inQ/l+O-06 +0.15 +0.24 
1U.O^_0.07 -0.16 -0.26 



Table III: Parameter point a at SPS la' and SPS 3 with la, 2(j, and 3a uncertainties (46 d. f.) from fits with 
the "LHC±ILC" setup. 



V. DISCUSSION AND CONCLUSIONS 



In this paper we have investigated the question how well one can measure i?-parity violating 
couplings using neutrino data and future collider data. For simplicity we have taken bilinear i?-parity 
violating parameters and fixed the J?-parity conserving ones by imposing mSUGRA boundary 
conditions. The latter choice is not important as the requirement of explaining correctly neutrino 
data only fixes ratios of i?-parity violating parameters over i2-parity conserving parameters. 

For the fits we have taken the current experimental accuracies on neutrino data and the expected 
accuracies on the measurements of edge variables at the LHC and mass measurements at the ILC as 
given by the corresponding colhder studies. In addition we have assumed that at the LHC the decay 
length of the lightest neutralino can be measured within 15 percent and that the branching ratios 
can be determined within twice the corresponding statistical uncertainties. Under this assumptions 
we find that the expected accuracies on the RPV parameters is of order one percent. However, all 
fits show that there will be a sign ambiguity as all observables considered are nearly the same under 
a sign change of the RPV parameters. This ambiguities might be resolved by detailed studies of 
the lepton and jet spectra of the individual decay channels of the lightest supersymmetric particle, 
in the case under study the lightest neutralino. 

The main source of the uncertainties are due to the uncertainties on the neutralino branching 
ratios. This implies that in a top-down approach, e. g. fitting the i?-parity conserving within a given 
high scale model such as mSUGRA, the main information is already obtained at the LHC even if 
the ILC measures the SUSY spectrum more precisely. However, improvements are expected if at 
at ILC the neutralino branching ratios can be measured significantly better than at the LHC as 
one would naively presume. In case one fits the i2-parity conserving parameters at the electroweak 
scale, then our findings hold if the masses of neutralinos, charginos, staus and sbottoms as well as 
the corresponding mixing matrices can be determined precisely. To keep things at the level shown 
(up to about a factor 1.5-2) one needs in case of the neutralino and charginos sectors precision in 
the percent-range, in stau and sbottom sector 10 % accuracies implying the need of an ILC and 
most likely also a multi-TeV e"^e~ collider such as CLIC. 
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